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Abstract 

We present the result for the three-loop singlet QCD corrections to the decay of a 
Higgs boson into two photons and improve the calculation for the non-singlet case. 
With the new result presented, the decay width T{II — t- 77) is completely known 
at 0{GFa'^aj,GFa^). 
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1 Introduction 



After the discovery of a new boson at about 126 GeV at the LHC pL|f2] it is 
imperative to study all decay modes of the Standard Model Higgs boson as 
precisely as possible. In this paper we are concerned with the decay of the 
Higgs boson into two photons which is the decay channel with the largest 
significance for the discovery of the Higgs boson at the LHC. 

The decay of a Higgs boson into two photons is mediated through either 
charged gauge bosons or top quarks at one-loop order. The decay rate can be 
cast into the form 

T{H^^^) = -^\A^ + A\', (1) 



Preprint submitted to Elsevier 



27 December 2012 



with the leading order values 
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where tw = Mfj/{4Myy) and r = Mfj/{Amf). The former contribution is 
larger than the latter (by a factor 4.5) and opposite in sign. Both contri- 
butions have been investigated in great detail in the literature. The two- 
loop QCD corrections to the decay have first been evaluated in the heavy- 
top limit in Refs. [3P|5] and later, keeping the full top-mass dependence, 
in Refs. |6]I71I8] . The two- loop electroweak corrections have been investigated 
in pilOm 11112] . Combining the two-loop QCD and electroweak corrections one 
observes a nearly complete cancellation between these two contributions for 
Mh = 126 GeV as discussed below. 



At next-to-next-to-leading order (NNLO) the non-singlet QCD contributions 
(cf Fig. [1] (c)-(e)) have been calculated in the heavy top hmit, including ad- 
ditional terms in an expansion in r = Mfj/iAmj), in Ref. [13]. At this order 
a new class of diagrams, so-called singlet diagrams (cf Fig. [1] (f)-(g)), con- 
tributes for the first time. They can be characterized by the property, that the 
external lines are coupled to different fermion loops. The contribution from 
this kind of diagrams has not been taken into account up to now and there is 
no formal argument, that it should be suppressed compared to the non-singlet 
contribution, when considering the same order in a^. In this letter we present 
the calculation of this last missing piece to obtain a complete NNLO QCD 
prediction for the decay of a Higgs boson into two photons. In addition, to 
improve the existing prediction and to check our setup, we recalculated the 
non-singlet contribution and added more terms in the expansion in r. 



This letter is organized as follows: In Section [2] we describe the necessary 
steps of the calculation and in Section [3] we present the full three-loop QCD 
corrections to the decay width. In Section H] we give the numerical result and 
our conclusions. 
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Fig. 1. Sample diagrams of the 1-loop (a), 2-loop (b), 3-loop non-singlet (c)-(e), and 
3-loop singlet (f)-(g) top quark induced contribution to — )• 77. 

2 Calculation 



We start from the decay amplitude of a Higgs boson into two photons 



A 




(3) 



Decomposing the Lorentz structure in the most general way the amplitude 
A^"^ can be written as 



(4) 



with scalar form factors A, B, C, D, and E, which can be calculated by 
applying the projectors 
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on the amplitude (e. g. A = P'^ A^y). Due to gauge invariance the contractions 
qinA^" and q2vA^^ of the amplitude with the momenta of the external photons 
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vanish. This imposes the relations B = —qiq2A and D = £" = on the form 
factors. In order to verify these relations explicitly we calculate A, . . . , E 
separately. The coefficient C, although non zero, does not contribute to the 
decay rate of the Higgs boson because of the transversahty of the photon wave 
functions. Therefore the Lorentz structure of the physical amplitude simplifies 
to 

= {qiq^-qiq2gnA. (6) 

To obtain the decay width, the combination A = A^ + At is inserted into 
eq. ([1]) for the W and the top quark loop induced contributions. 

Since the calculation of the decay amplitude including the full dependence 
on Higgs-boson and top-quark mass is not possible at the moment, only an 
expansion in r = Mfj/{4mf) can be obtained. The expansion of the amplitude 
in T can be performed following two conceptually different strategies. The first 
option, which is only applicable for the calculation of a few terms in the ex- 
pansion, relies on the direct expansion of each Feynman diagram in r and the 
subsequent reduction of the appearing integrals to master integrals. Since this 
method becomes very time-consuming when calculating higher terms in the 
expansion, a second approach has been devised. In this alternative approach, 
the full mass dependence is kept as long as possible, meaning that the reduc- 
tion to master integrals is performed keeping the full mass dependence and 
only the master integrals are calculated in an expansion in r. In the following 
we sketch the steps necessary to perform the calculation and explain how the 
results for the needed master integrals in an expansion in r can be obtained. 

The Feynman diagrams are generated with qgraf [T3] and processed by a 
Mathematica program to map the integrals on 4 singlet and 15 non-singlet 
topologies. A sample of the contributing diagrams is depicted in Fig. [TJ All 
appearing integrals are reduced to master integrals using Crusher [L5\. To 
expand the master integrals Mj(r) in r = Mfj/lAmf) a series ansatz 

M^ir) = Y: + Ml'^i-rr + Ml^i^rY^y (7) 

k 

is inserted into the differential equation [16j 

where ID applied to Mj(r) returns the mass dimension of Mj. The differen- 
tial equation follows from the fact that Mj(r) is a homogeneous function in 
its dimensionful parameters Mjj and m^. Quarks except for the top quark 
are treated as massless. Using the series ansatz results in a system of alge- 
braic equations which are systematically solved for the coefficients Mf^^ by a 
Mathematica program which uses Fermat [T7] to simplify rational functions. 
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Fig. 2. Integrals which serve as boundary conditions for the solutions of the differ- 
ential equations for the vertex master integrals. Dashed internal lines are massless, 
solid lines carry the mass m^. External lines on the right side of the diagrams (a)-(e) 
are massless on-shell, the external line on the left side is on-shell with mass Mh- 
These diagrams only contribute to the singlet case. 

By this procedure all coefficients Mj^'* are expressed as linear combinations of 
the integrals which are depicted in Fig. [2] and serve as boundary conditions. 

(k) (k) 

Coefficients M>1' and Mf^' of non-integer powers of r arise only in the singlet 
contribution due to the massless cuts of the singlet diagrams [18] . 



3 Results 



The top quark loop induced amplitude At is expressed as a perturbative se- 



ries in the strong coupling constant with the one-, two-, and three-loop 

2 



contributions A^ , A[ , and 



A = A 



71 V vr / 



(9) 



The three-loop contribution is, furthermore, split into the non-singlet part 

(2) (2) 

AIq, the singlet part with two top-quark loops A\J , and the singlet part with 

(2) 

a top-quark and light-quark loop AIJ 



4^) = Ag + 4J + Wr^EQ?)^S (10) 



Qt and Qg are the electromagnetic charges of the top quark and the light 
quarks (g G {u, d, s, c, b}). 

The full results with the expansion in r up to r^'' including generic SU{N) 
colour factors and renormalisation scale dependence are available onlineQ For 
the presentation of the results in this section we truncate the expansion after 
and insert the SU{3) colour factors explicitly. The top-quark mass is renor- 
malized in the MS scheme. We use the abbreviations = log fi^/mf in the 



^ http : //www-ttp .particle .uni-karlsruhe . de/Progdata/ttpl2/ttp 12-046/ 
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non-singlet and = log(— 4r) in the singlet contribution. For completeness 
also the and al contributions are given as an expansion in r. 
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.(2) 13 2 ^ 1 ^ 

At a = + -Cs + -Lr 

t,q 12 3^^ 6 
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The results for q agree with the results presented in Ref. [13]. In Tab. [T] 
we give the results for the first 20 coefficients of the series in r in numerical 
form. For the singlet contribution we show the results for the constant and the 
logarithmic part proportional to Lj. separately. To illustrate the convergence of 
the series we show all 20 terms of the expansion of the three-loop contribution 
in graphical form in Fig. |3l Around r = 0.14, corresponding to the Higgs- 
boson mass favoured by the LHC experiments, the first four (five) terms in 
the expansion are sufficient to obtain an accurate result with a relative error of 
10~^ wrt. the total 3-loop amplitude for the singlet (non-singlet) contribution. 
At larger values of r more terms of the expansion have to be taken into account, 
e.g. at r = 0.5 nine (sixteen) terms are needed to obtain similar accuracy. 

Collecting all available information the amplitude can be cast in the form 

^LO ^NLO-EW ^NLO-QCD J^nnLO 

where denotes the electroweak corrections to the W and top-quark in- 
duced processes combined. The partial decay width is then given by 

^H^-yy = 7^;— l^LO — (2v4lo^nlo-ew) H (2Alo^nlo-qcd) 

o47r ^ vr tt 

+ (^^) (2ALoRe(ANNLo) + ^NLo)) • (17) 

For a Higgs boson with a mass of Mh = 126 GeV (r = 0.14) this evaluates 
numerically to 

Th^^^ = (9.398 ■ 10"^ - 1.48 ■ 10"^ + 1.68 ■ 10"^ + 7.93 • 10"^) GeV (18) 
= 9.425 ■ 10"^ GeV, 
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T 

Fig. 3. Numerical results for the non-singlet and singlet contribution at three-loop 
order. In the plot the first 20 terms in the expansion in r are shown. The singlet plot 
shows the sum a[^^ + Sg^t Qg/Qt^Lq- ^'^^ the plot of the non-singlet contribution 

(2) 

A^Q we set the renormalization scale /i = Mh- The vertical line at r = 0.14 marks 
the value of r for Mh = 126 GeV. The convergence of the expansion for r — )• 1 can 
be improved by using the on-shell scheme for the top-quark mass. 

where we used mt(Mjy) = 166 GeVH, a,(Mj^)/7r = 0.0358, Gp = 1.16637- 
10^^ GeV^^, and a = a{0) = 1/137 as input parameters. The value for the 
two- loop electroweak correction was taken from Ref. [TU]. It has to be noted, 
that there is a partial cancellation between the two-loop QCD and electroweak 
corrections. To assess the influence of the singlet diagrams the next-to-next- 
to-leading order term can be further decompsed as 

r_f/^77|NNLO = (7.5 ■ 10 ^°|nlo2 + 1-73 ■ 10 ^|non-singlctt + 5.45-10 ^ | singlet) GcV , 

which shows, that the singlet diagrams are a factor of three larger than the 
non-singlet ones and thus the most important three-loop contribution. 



4 Conclusion 

We presented new results for the singlet diagrams which contribute to the 
decay of a Higgs boson into two photons at next-to-next-to-leading order, 

^ The MS mass mt{M}j) was obtained from the on-shell mass Mt = 172.64 GeV 
using the MS-on-shell relation at 3-loop accuracy |19|20|21] . 
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Tabic 1 

Numerical values of the first 20 terms in the expansion in r for the three-loop non- 

(2) (2) (2) 

singlet A), q and singlet contribution Al J and A\^q . For the non-singlet part we set 

' (2) 
IJ, = uit-ln case of the singlet contribution the results for the constant part {Ai [0]) 

and the logarithmic part are shown separately. 

which up to now have not been considered. The corrections to the decay rate 
due to singlet diagrams are about a factor of three larger than the non-singlet 
ones. An improved prediction for the non-singlet contribution reduces the 
error of the three-loop contribution, which can be neglected for a Higgs-boson 
mass of 126 GeV. The total partial decay width is T//^^^ = Tlq + Tnlo-ew + 
rNLO-QCD + Tnnlo = (9.398 ■ 10-6 - 1.48 • 10"^ + 1.68 ■ 10-^ + 7.93 • IQ-^) GeV = 
9.425 • 10-6 GeV. 
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